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Abstract 

vC . 

Q\ , We use the wounded nucleon model to study the forward-backward mul- 

tiplicity fluctuations measured by the PHOBOS Collaboration in Au + Au 
collisions at ^/s NN = 200 GeV. The enhancement of forward-backward fluctu- 
ations in Au+Au collisions with respect to the elementary p+p interactions is 
in this model explained by the asymmetric shape of the pseudorapidity density 
of produced particles from a single wounded nucleon and the fluctuations of 
the number of wounded nucleons in the colliding nuclei. The wounded nucleon 
model describes these experimental data better than the HIJING, AMPT or 



X 

UrQMD models do. 
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1 Introduction 

The studies of the fluctuations of the charged particles multiplicity in heavy nuclei 
collisions provide information on the particle production and the properties of the 
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matter created in these collisions pQ. The multiplicity measured in the symmetric 
pseudorapidity intervals in the forward and backward hemisphere appears to be a 
suitable observable for such analysis. Recently, the direct forward-backward corre- 
lations were studied by the STAR Collaboration |2j, while the difference between 
forward-backward multiplicities was used by the PHOBOS Collaboration [3j. 

The PHOBOS Collaboration measured the variance Oq of the forward-backward 
asymmetry variable C = (N F — N B )/\ / N F + N B 



where N B and Np are the event-by-event multiplicities in the backward B and 
forward F pseudorapidity intervals, respectively. The measurement was performed 
for many different symmetric (with respect to rj = 0) pseudorapidity bins in the 
range \rj\ < 3. It was found that: 

(i) for rj bins with fixed width but changing position in pseudorapidity o 2 c in- 
creases with increasing distance between bin centers 

(ii) for T) bins centered at the same position [r\ = ±2) a 2 c increases with increasing 
width of the bins 

(iii) Oq is significantly larger for peripheral (centrality 40 — 60%) collisions than 
for central (centrality — 20%) collisions. 

In this paper we use the wounded nucleon model [3] to calculate the values of 
Oq for the the same pseudorapidity bins as in the PHOBOS experiment. First we 
perform calculations for the simple case of elementary p + p interactions and then 
we extend them to Au + Au collisions and different centralities. We show that the 
wounded nucleon model reproduces reasonably well the experimental data, especially 
the centrality dependence of the forward-backward fluctuations. 

In the next section we define the model in detail. In Section 3 we describe 
our Monte Carlo calculations and present the approximate analytical formula for 
<Jq derived in the framework of the wounded nucleon model in the Appendix. In 
Section 4 our results are compared with the PHOBOS data and also predictions 
of other models of nucleus-nucleus collisions are discussed. We end our paper with 
conclusions. 

2 Wounded nucleon model 

The wounded nucleon model [3J E] assumes that after the collision of two nuclei, the 
secondary particles are produced in the process of independent fragmentation^ of 

is similar to the assumption of independent hadronization of strings in the dual parton model 




(1) 



6j. 
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Figure 1: The wounded nucleon fragmentations function p{rj): pseudorapidity den- 
sity of produced particles from the forward-moving wounded nucleon at ^fs = 200 
GeV cm. energy. 

these nucleons (called "wounded nucleons"), which underwent at least one inelastic 
collision. The contribution from a wounded nucleon does not depend on the number 
of such collisions^!, thus the secondary particles are produced according to a universal 
fragmentation function, p{rj). Number of produced particles may vary following a 
multiplicity distribution. 

The shape of the fragmentation function was extracted [SI H\ from the PHOBOS 
data on d + Au collisions [8] at ^/i = 200 GeV and is presented in Fig. [1] for a 
wounded nucleon from the forward moving nuclejE The fragmentation function for 
nucleons from the opposite nuclei is a reflection with respect to rj = 0. 

In this model the single particle density in elementary p + p collisions is de- 
scribed as pijf) + p{—rj). In the case of nucleus-nucleus collision it becomes a sum 
of these functions multiplied by the number of wounded nucleons in the appropriate 
nuclei [S]: (wf) p(t)) + (wb) p(—r)), where (wp) and (wb) are the mean values of 
wounded nucleons in the forward and backward moving nucleus, respectively. In 
consequence at midrapidity the particle density in A + A equals (w) p(0) (where 
(w) = (wf) + (wb)) and is thus (w) /2 times larger than in elementary collisions. 
This prediction does not agree with the experimental data, however, as the parti- 
cle density dN/dr]\ v=0 is larger than (w) p(0) [TO]. This can be easily corrected by 
introduction of a centrality dependent factor 7 which appropriately increases the 
multiplicity of particles generated in the fragmentation of a wounded nucleon that 

2 In fact in our calculations we will modify this assumption. We will come back to this point 
later. 

3 An analogous wounded nucleon fragmentation function was constructed in Ref. [S]. 
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underwent many inelastic collision^. 

The necessity of the above correction results from the composite structure of 
the nucleons. It is accounted for in the wounded quark-diquark model [7j, where 
the fragmenting objects are quarks and diquarks, and which reproduces the particle 
density measured in A + A collisions. The fragmentation of a wounded nucleon 
is equal in this approach to a sum of fragmentation functions of the quark and 
diquark weighted by the probability of " wounding" each of them. The scaling factor 
7 introduced by us reflects the changes of these probabilities for different centralities 
of the collisions. The wounded nucleon model with our correction is equivalent to 
the wounded quark-diquark model in the kinematical area, where the fragmentation 
of spectators (unwounded constituents) is negligible [TJ. This is fully justified in our 
case, as we are analyzing data restricted to \rj\ < 3. 

The fragmentation function p(rj) represents the mean yield from a wounded nu- 
cleon, whereas for calculation of fluctuations the changes occurring event-by-event 
have to be included. We assume that the shape of pirf) is unchanged, even if the 
integral of p(rj) follows a multiplicity distribution. Further we also assume that the 
multiplicity distributions in p+p collisions in the combined pseudorapidity bin B + F 
can be parameterized by negative binomial (NB) distribution [11] 



T(n + k) m\ n / n\- n - k 

p <"'"'*> = r(„ + i)r W (t) + • (2) 

where n is the average multiplicity in the combined interval B + F, 1/k measures 
deviation from Poisson distribution and T(n) is the gamma function. In our calcula- 
tions we need to know the multiplicity distribution from a single wounded nucleon. 
Assuming that soft particle production in p + p collisions can be described by inde- 
pendent contributions from each nucleon [12J, it is straightforward to show that the 
multiplicity distribution from a single wounded nucleon follows also NB distribution 
(j2J) but with n and k substituted by n/2 and k/2, respectively^. 

To perform the necessary calculations we also need to know the probability that 
a particle originating from the forward moving wounded nucleon goes to F interval, 
under the condition it is found either in B or F. It can be easily calculated as 

p f F p(v) d v (3) 

f B p{rj) drj + f F p{rj) dr]' 

where p (rj) is the pseudorapidity density of produced particles shown in Fig. [U 
In consequence the probability that a particle originating from the forward moving 



4 Contrary to the wounded nucleon model assumption multiplicity from a wounded nucleon 
slightly depends on the number of collisions it underwent. 

5 In general the multiplicity distribution from w independent wounded nucleons is given by NB 
distribution @ with h — > wn/2 and k —> wk/2, where n and k are measured in p + p collisions. 
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wounded nucleon goes to B interval equals I— p. For the backward moving wounded 
nucleon p has the same meaning when the intervals F and B are exchanged. The 
n parameter from NB distribution can be also calculated using the fragmentation 
function 



For the nucleus-nucleus collisions in our modified version of the wounded nucleon 
model the multiplicity of particles from a single wounded nucleon is rescaled by 7 
factor. It reflects the fact that a wounded nucleon which underwent many inelastic 
collisions produces more particles than a wounded proton in p + p collisions. Ac- 
cordingly, in Eq. (T5]) we substitute n not by n/2 but by 771/2. Less obvious is the 
necessity of modification of k parameter, which we decided to multiply also by the 
factor 7. This is justified if we assume that each nucleon is composed of two con- 
stituents, a quark and a diquark, which populate particles independently. However, 
we show also results for unmodified k and we consider the difference between them 
as an additional systematic uncertainty of the wounded nucleon model. 

3 Calculations 

In the present section we describe the MC calculation of <Jq defined in Eq. ([1]) 
in the modified wounded nucleon model. In addition we present also approximate 
analytical formulae (with n and k rescaled be the factor 7) derived in the Appendix. 
We performed our calculations in three steps described below. 

(A) In order to study effects of simple superposition of nucleon-nucleon collisions 
we assume that in both nuclei we have fixed number of wounded nucleons wp = wb- 
Each wounded nucleon adds it's contribution generated independently according to 
the NB distribution (Eq. [2]) with parameters n and k substituted by 7n/2 and jk/2 
respectively, where n and k are measured in p + p collisions. The generated particles 
are then assigned randomly, with probability p and 1 — p, to the intervals F and 
B (for the forward moving wounded nucleon; to B and F for the backward moving 
nucleon) . 

Our main observation is that this contribution very quickly saturates with in- 
creasing wp. For wp > 10 it practically does not depend on wp. In the Appendix 
we present the approximate analytical derivation of the asymptotic formula for large 
Wf, which is: 



This term depends only on the ratio n/k, thus it is insensitive to our modification 
of the wounded nucleon model, i.e. the simultaneous substitution n i— > 777, and 




(4) 




(5) 



5 



k i — > jk. Except the most peripheral events 0^| A is centrality independent. 

(B) Contribution solely from the fluctuations in the number of wounded nucleons 
can be obtained when we assume that each wounded nucleon fragments always into 
the same, mean number of particles. Then each forward moving wounded nucleon 
populates 'ynp/2 particles into the forward F interval and 7^(1 — p)/2 particles into 
the backward B interval (and analogously for the backward moving nucleons). In 
consequence we obtain 

N F = \pn [w F p + w B (l- p)} , 

n b = ]pn [w B p + w F (l -p)} , (6) 

where w F and w B are the numbers of wounded nucleons in the forward and backward 
moving nuclei, respectively. The above formulae allow to substitute N F and N B in 
the definition of (Eq. [[]): 

«| n = i7"(2p-W^^y (7) 
2 \ w F + w B 



This term explicitly depends on w F and w B and thus is sensitive to the centrality 
of the events sample. Moreover this term is sensitive to the 7 factor which was 
introduced to correct the wounded nucleon model. The distribution of w F and w B 
cannot be obtained in a model independent way. In order to match the experimental 
centrality definition we are using the events from the HIJING generator, for which 
cuts compatible with these for the PHOBOS experimental data were performed [T3] . 
However, similar results can be obtained if w F and w B are obtained from the Glauber 
type MC generator using centrality cuts on impact parameter or on % + w^. 

(C) The full calculations combine (A) and (B) together and include both the 
multiplicity distribution and the fluctuations of the number of wounded nucleons. 
The w F and w B are taken from generated HIJING events. Next each wounded 
nucleon populates particles according to the distribution ([2]) in a procedure similar 
to that from (A). 

Results of full simulations are practically identical with the sum of two previously 
described terms (A) and (B). This is confirmed by the analytical calculations in the 
Appendix as in the first approximation (neglecting terms of the order 1/ (wf)) we 
obtain 



1 + n(2p - 1 



,2 



11/ (wp — w B 
>7 



k 2 \ w F + w B 



(8) 



6 We checked it with our own MC generator with the Gaussian nucleon-nucleon interaction 
profile. 
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Table 1: The parameters used in the calculations of wounded nucleon model predic- 
tions presented in Figs. [2] and [31 We estimate the uncertainty of the values of the 
parameters to be 7%, 20% and 10% for An/n, Ak/k and Ap/p, respectively. 



According to this formula, the values of Jq > 1 can be obtained in the wounded 
nucleon model only when the wounded nucleon fragmentation function is asymmetric 
with respect to rj = (leading to p ^ 0.5). The effects increasing a^, already 
contained in the multiplicity distribution measured in elementary p + p collisions 
are enhanced by superposition of nucleon-nucleon collisions. The fluctuations of the 
number of wounded nucleons, that cause the difference wb — Wf to be non-zero, 
further increase the a 2 c . 

This closes theoretical discussion of the problem. 



4 Results and discussion 

In the present section we compare results for the variance a 2 c of the forward- 
backward asymmetry variable C = (N F — N B ) / y/Np + Nb obtained in our wounded 
nucleon model MC simulations with the recently published PHOBOS data [3j. The 
measurement was performed for Au + Au collisions at ^/s = 200 GeV in the pseu- 
dorapidity range of \r]\ < 3 for various symmetric forward and backward intervals. 

The model calculations were performed for the same intervals as in the experi- 
ment and for the experimentally studied centrality classes. The parameters p and 
n were obtained from the wounded nucleon fragmentation function shown in Fig. [1] 
using Eqs. (J3j) and (HI) and are listed in Tab. [TJ The 7 factor responsible for the 
correction of the wounded nucleon model was estimated by comparing (w)p(0) and 
the PHOBOS data on dN/dr] at midrapidity [10], and equals 7 = 1.35 ± 0.15 and 
1.6 ± 0.1 for peripheral and central collisions, respectively. 

The dependence of a 2 c on the bin width A77 presented in Fig. [2] includes points 
with the largest values of a^, measured with the smallest relative errors, and is thus 
more challenging for the models. Here the center of each bin is fixed at rjp = 2 and 
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Figure 2: The variance of the forward-backward asymmetry variable C = (Np — 
Ng) I \J Np + Nb as a function of the width of the forward and backward intervals 
for most central collisions (left) and the peripheral collisions (right). The position 
of the center of each bin is fixed at rjp — 2 and rjs = —2. The PHOBOS data 
points (dots) are compared with the prediction of the modified wounded nucleon 
model (continuous and dashed lines, for the scaled and unmodified k parameter 
respectively). The grey bands reflect the systematic error of our calculations, due 
to the uncertainty of the parameters (see the text near the end of Section H]). 

i]b = —2. Our modified wounded nucleon model correctly reproduces the general 
trends observed in the data, especially the centrality dependence. The predictions 
obtained for the case with scaling of k are slightly below the experimental points, 
but in most cases are consistent within the systematic errors. Interestingly, in the 
case of unmodified k, which we consider less probable, the agreement is much better. 

In Fig. [3] the values of a 2 c obtained for fixed width of each interval (At] = 0.5) 
and changing the distance between bin centers are shown. The wounded nucleon 
model is in this case in worse agreement with the data, there are discrepancies 
for 0.5 < 7] < 2. However, for narrow rj bins the values of a 2 c are closer to 1, the 
systematic errors are relatively larger and also the influence of effects not included in 
the model may be thus stronger. The NB distribution which is supposed to contain 
effectively all correlations present in the elementary collisions may be not sufficient 
to fully represent the short range correlations (resonance decays, clusters). 

In Fig. H] we show different components of a 2 c present in the wounded nucleons 
model (with scaled k parameter), as discussed in the previous section. For elemen- 
tary p + p collisions the value of Oq for the widest interval (Arj = 2) reaches 1.34, 
the superposition of nucleon-nucleon collisions increases it to about 1.5. There is 
almost no difference between wp = 30 and wp = 140, which correspond to the mean 
values for the peripheral and central collision analysed, as the superposition effects 



S 




Figure 3: The same as in Fig. [2] but now the width of each bin is fixed and equals 
At] = 0.5, changing 77 is the position of the center of the forward bin (the backward 
one is symmetric around 77 = 0). 

fast saturate. Much stronger increase is due to contribution from the fluctuations 
of the number of wounded nucleons, they add to the previous numbers directly and 
give 1.94 and 2.41 for central and peripheral collisions respectively. 

Let us notice that all presented results^ can be obtained from our analytical 
approximation (jSJ), which agrees with the presented MC results with the accuracy 
better that 3%. In the case of analytical calculations we are using the values of 
^ ^wp+wb ) ex t r acted from samples of peripheral and central HIJING events selected 
according to the experimental centrality cuts, which are 0.38 ± 0.05 and 0.90 ± 0.09 
for centrality — 20% and 40 — 60%, respectively [T3] . 

Finally, let us discus the uncertainty of our approach, represented by the grey 
bands in Figs. [2] and [31 The systematic errors of model calculations can be es- 
timated using the analytical formula (jSJ). We calculate the error of caused by 
the uncertainty of each parameter separately and add them in quadrature. In the 
case of k parameter the systematic error contains only the error of k (20%), the 
uncertainty of scaling k for A + A collisions is represented by two extreme cases, 
for which the systematic errors are calculated separately. All contributions due to 
errors of parameters are of the same order, but the uncertainty of the probability p 
has always the largest impact. The systematic error of predicted a 2 c calculated this 
way is approximately equal to 0.17 (o^ — 1). 

The data analyzed in this paper were previously compared with expectations 
from several models of nucleus- nucleus collisions: HIJING [15], AMPT [16] in Ref. 
[3] and UrQMD in Ref. [H]. As we can see in Fig. [5] the differences between them 



7 Except the predicted values of u 2 c for elementary p + p collisions presented in Fig. [3J 
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Figure 4: Wounded nucleon model predictions (with scaled k parameter) for 
values obtained for the pseudorapidity bins centered at rj = 2 with changing width 
Arj. The model calculations for central 0—20% ((%) = 140) and peripheral 40—60% 
{(wp) = 30) Au + Au collisions with all effects included (full and open circles) are 
compared with results for collisions with exactly 140 and 30 wounded nucleons in 
each nuclei (full and open triangles), the fluctuations in the number of wounded 
nucleons is responsible for the difference between corresponding sets of points. In 
addition the predictions for p + p collisions are also shown (stars). 

are large, but none of these models correctly describes all features of the data listed 
in the Introduction. HIJING gives approximately correct values of <j 2 c for central 
collisions, but they do not grow for peripheral collisions. AMPT features expected 
centrality dependence, but always significantly underestimates o 2 c values. UrQMD 
calculations agree with the experimental results for peripheral events, but this model 
predicts a slight increase of afj for central collisions in place of an observed significant 
drop. The predictions of the wounded nucleon model are closest to the experimental 
data - especially for the version with k parameter not scaled. 

5 Conclusions 

Our conclusions can be formulated as follows. 

(i) We show that the PHOBOS Collaboration data [3] on the variance Oq of 
the forward-backward asymmetry variable C = (N F — N B ) / \J Np + Nb measured in 
Au + Au collisions at \/s = 200 GeV can be reasonably described in the framework 
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Figure 5: Comparison of the experimental data with predictions from several models 
for (j^ dependence on the width A?? of the pseudorapidity bin. Results for wounded 
nucleon model are shown for two cases, k parameter unsealed and scaled, as dashed 
and continuous lines respectively. UrQMD [14], HIJING and AMPT [3] predictions 
are represented by open symbols described in the plot. 



of the wounded nucleon model. The wounded nucleon model predictions in most 
cases are consistent within the errors with the data and slightly lower values of a\, 
may be due to short range correlations, which are only partially accounted for in the 
model. Anyway, we reproduce both the pseudorapidity and centrality dependence 
of Oq in contrast to much more advanced HIJING, AMPT or UrQMD models. 

(ii) The key ingredients of our approach are: the assumption of independent 
hadronization of each wounded nucleon; the wounded nucleon multiplicity distribu- 
tion described by the negative binomial distribution (T5]) and; the wounded nucleon 
pseudorapidity fragmentation function, shown in Fig. [TJ Moreover, we correct the 
wounded nucleon model by taking into account the dependence of the multiplicity 
from a single wounded nucleon on the number of collisions. 

(iii) We also derive approximate analytical formula for cr^, (jSj). It can be observed 
that o'q — 1 is proportional to (1 — 2p) 2 where p determines the partition of each 
particle between forward and backward intervals. Thus the PHOBOS data can be 
reasonably described only if the wounded nucleon fragmentation function is indeed 
asymmetric in pseudorapidity, i.e., p ^ 0.5. It is also interesting to note that the 
centrality dependence of a% is fully determined by the simple factor ( ( w b-*>f) \ 

Acknowledgements 

We thank the PHOBOS Collaboration for the access to the simulation data which 
enabled to calculate the model predictions most compatible with the experimental 
conditions. Discussions with A. Bialas are highly appreciated. This investigation 



11 



was supported in part by the Polish Ministry of Science and Higher Education, 
grant No. N202 125437, N202 034 32/0918, N202 282234 (2008-2010) and by the 
Director, Office of Energy Research, Office of High Energy and Nuclear Physics, 
Divisions of Nuclear Physics, of the U.S. Department of Energy under Contract No. 
DE-AC02-05CH11231. A.B. acknowledges support from the Foundation for Polish 
Science (KOLUMB programme). 



A Appendix: Calculation of 

In this Appendix we present the approximate analytical calculation of a 2 c . Let us 
introduce the following notation 

N B = (N B )(l + e B ), 

N F = (N F )(l + e F ). (9) 

Although generally e B and e F may have any value from —1 to infinity, if we exclude 
the most peripheral events, the values of \e B \ and \e F \ are much smaller than 1. In 
such a case it is possible to expand the denominator of cx^. in the powers of e B and 
e F (here we consider collisions of identical nucleus and symmetric B and F bins, i.e. 

(N B ) = (N F )) 

1 1 



N B + N F (N B ) \2 + (e B + e F 



1 fl e B + e F | je B + e F ) 2 | . ,. |( ^ 



(N B ) \2 4 

Taking the first terrrj§ of this expansion into account we obtain much simpler formula 
to handle 

(Nb - N F f \ _ (N B (N B - 1)) - (N B N F ) M 
N B + N F / ~ (N B ) 

Indeed, now we can apply exactly the same method as in Ref. [17J, where the 
forward-backward correlation coefficient in the wounded nucleon model was calcu- 
lated. Let us construct the generating function 

H(z B ,z F ) = P(N B ,N F )z^ B z^, (12) 

N B ,N F 

where P(N B , N F ) is the probability to find particles in the backward B interval 
and N F particles in the forward F one. Of course P(N B , N F ) depends on the 



5 Calculation of the higher orders is straightforward. 
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centrality of the collision. In Ref. [17] we have shown that in the framework of the 
wounded nucleon model, defined in Section 2, we obtain 

y| > -kw B /2 

X 



H(z B ,z F ) = J2 W(w B ,w F ){l + ^\p(l-z B )+lf(l-z F )]} 

x[l + -\p'(l-z B )+p(l-z F )}j , (13) 



where W (w B ,w F ) is the probability distribution of the numbers of wounded nucle- 
ons in the backward and forward moving nucleus, respectively, p' = 1 — p (B and 
F are symmetric around r\ — 0) where p is defined in Eq. Q. n and k come from 
the NB fits (J2]) to the p + p multiplicity data in the combined bin B + F . Now 
the averages present in (fTTT) can be directly calculated using the generating function 

am 



I at \ dH(z B ,z F ) 



(N B (N B -1)) 
(N B N F ) 



dz B 




d 2 H(z B , 


z F ) 


dz% 




d 2 H(z B , 


Zp) 



dz B dz F 



zb=zf=1 
zb=zf=^ 

(14) 

zb=zf=1 



Performing appropriate differentiations, multiplying n and k by the 7 factor and 
taking the following relation into account 

E w {WB , mF) ^pf = (J^pfl „ /(sz^^n (15) 

^ 2 (w B ) 2 (wb) \ w B + w F / 

Wg,WF 1 ' 1 ' v ' 

we obtain our final result (151). 
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